10 PRINT "USE OF SIMPSON'S RULE TO INTEGRATE G(Y) FROM 1 TO 2": PRINT
20 DEFINT I-J, N: DEF FNG (YY) = 3.5 * YY + LOG(YY) + 40! * SIN(1.57 * YY)
30 N = 20: Y0 = 1!: Y1 = 2!: H = (Y1 - Y0) / N
40 V = 0!: N1 = N - 1: FOR J = 1 TO N1 STEP 2: I = J - 1
50 V = V + H / 3 * ((FNG(Y0 + I * H)) + FNG(Y0 + (I + 2) * H) + 4! * FNG(Y0 + (I + 1) * H))
60 Y = Y0 + (I + 2) * H: PRINT "  Y="; Y; "  V="; V; : PRINT : NEXT J: PRINT